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I.  INTRODUCTION  AND  SUMMARY 

As  a  homing  missile  closes  on  its  target,  a  larger  portion  of  Its  seeker 
beam  Is  subtended  by  the  finite  angular  extent  of  the  target.  The  radar  echoes 
from  the  target  no  longer  appear  to  arrive  from  a  single  direction.  The 
direction  to  the  target  is  determined  from  the  measured  amplitudes  and  phases 
of  the  received  difference  and  sum  channel  signals  as  j^-jcos  (<pD  -  <pg)» 

Mutual  interference  between  the  echoes  from  individual  scatters  causes  the 
phase  front  across  the  antenna  (and  hence  the  apparent  direction  to  the  target) 
to  oscillate  about  its  mean.  For  a  point  target  the  phasors  D  and  S  are 
collnear.  When  phase  front  distortion  occurs,  the  quantity  ^  sin  (<J>D  -  <Pg) 
is  no  longer  zero  and  the  sum  channel  amplitude  undergoes  a  fade.  The  apparent 
glint  in  the  resulting  target  direction  can  be  reduced  by  designing  a  signal 
processor  which  senses  the  occurrence  of  this  interference  phenomenon  and 
rejects  the  consequent  erroneous  tracking  data. 

This  paper  presents  the  analysis  of  such  a  processor  and  shows  the  glint 
reduction  achieved  by  applying  a  combination  of  thresholds  to  |s|  and 
j-^j  Isin  (<J>D  -  (j>s)  | .  Curves  are  presented  which  show  the  fraction  of  data 
rejected  by  the  thresholding  procedure  as  a  function  of  the  two  applied 
thresholds.  Two  measures  of  glint  reduction  are  evaluated;  curves  are  pre¬ 
sented  for  each  of  these  measures  which  show  the  glint  reduction  obtained 
versus  the  fraction  of  data  lost.  For  a  specified  fraction  of  rejected  data, 
maximum  glint  reduction  is  obtained  by  thresholding  on  |s|  alone. 

II.  ANALYSIS 

The  homing  missile  employs  a  phase  comparison  monopulse  seeker;  each 
subaperture  of  the  antenna  is  assumed  to  possess  an  isotropic  gain  pattern. 

We  consider  the  glint  effects  in  a  single  angle  dimension  only  (Figure  1). 


If  a  monochromatic  wave  is  reflected  from  a  point  target  at  range  R»d 
and  off-boresight  angle  ip  then  the  RF  signals  received  at  -  +  y  respectively 

are  given  by  . 
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Figure  2.  Monopulse  error  characteristic. 
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where  k'  contains  the  target  amplitude  and  target  phase  referenced  to  the 
transmitted  signal  from  the  Illuminating  radar. 

Magic  T's  at  the  antenna  combine  these  signals  so  that  the  receiver  Input 
consists  of  the  sum  and  difference  signals. 

S  A  V  (-|)  +  V(+  |)  (3) 

t 

D  A  +  1  [v(-  |)  -  V(+  |)]  (4) 

Since  R»d,  the  complex  envelopes  of  the  IF  signals  In  the  sum  and  difference 
channels  of  the  receiver  may  be  written  as 

S  *  A  cos  y  A_  u  +  iv  (5) 

D  -  A  sin  Y  jA  X  +  iy  (6) 

where  the  complex  amplitude  A  Includes  the  effect  of  range  attenuation,  target 
RF  phase,  target  cross-section,  etc. ,  and  we  have  denoted  by 

Y<*>  ■  ysln  i(/  (7) 

the  RF  phase  shift  for  a  target  at  $  referred  to  the  RF  phase  which  would  be 
received  at  the  antenna  center.  (For  an  active  missile  y  is  replaced  through¬ 
out  by  2  y) •  For  a  single  target  the  ratio  D/S  is  real  and  is  Independent  of 
the  target  properties,  A.  Thus  by  forming  D/S  the  receiver  obtains  tan  Y  from 
which  the  off-boreslght  angle  is  uniquely  determined  (Figure  2). 

Although  for  a  point  target 

f-f-*(reaD 


(8) 


receivers  are  designed  to  produce  the  output 


cos  (<J>D  -  <f>g)  -  +  ^  -  Re  (D/S) 

u  +  v 


(9) 


This  has  the  virtue  of  not  becoming  Indeterminate  when  either  the  phase 

of  A  is  0  or  it  (y-v-0)  or  the  phase  of  A  is  —  or  ~  (x*u-0) ;  for  a  point  target 

Re  D/S  gives  the  same  monopulse  error  characteristic  as  does  D/S. 

Consider  now  signals  returned  from  an  extended  target  which  occupies  the 

6  4 

segment  of  the  beam  from  to  The  complex  envelopes  in  the  sum 

and  difference  channels  of  the  receiver  may  be  written  as 


S 


D 
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2 

cos  y(ij>)  drjj  +  ng  u  +  iv 
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2 

sin  y(iJ0  d^  +  nD  ^  x  +  iy 


(10) 


(11) 


where  ng  and  are  independent  receiver  noises  in  the  two  channels. 

Denote  statistical  averages  over  the  noise,  che  amplitudes  of  the 
individual  scatterers  and  the  phases  of  the  individual  scatterers  by  an 
overbar.  Then  assuming  that  the  phases  of  the  Individual  scatterers  are 
Independent  and  uniformly  distributed. 


S  -  u  +  iv  -  0 

(12) 

D  =  x  +  iy  -  0 

(13) 

Since  x  and  u  involve  only  the  cosines  of  the  target  phases  while  y  and  v 
involve  only  the  sines  there  is  zero  correlation  between  any  in-phase  and  any 
quadrature  signals 

xy-uv-xv-yu-0  (14) 
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Further,  for  either  of  the  channels  the  average  in-phase  and  quadrature  powers 
are  equal 
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(15) 
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sin2yWd\f>  +  ~  [n^  A  a2 


(16) 


Also,  the  average  cross-channel  powers  are  equal 


ro+  7 


XU 


yv 


I  AW  | 2  sin 
i|>  -6 


yW  cos  yW  W  4  Pab 


(17) 


2  2 

Thus  2b  ,  2a  and  2  pab  denote  respectively  the  average  sum  channel, 
difference  channel  and  cross-channel  powers. 

The  extended  target  is  assumed  to  consist  of  a  large  number  of  mutually 
interfering  scatterers  whose  returns  generate  a  Gaussian  process. 

Thus,  (12)  .  (17)  imply  that  the  in-phase  pairs  x,  u  and  the 

quadrature  pairs  y,  v  obey  independent,  bivariate,  zero-mean  Gaussian  distri¬ 
butions.  Because  of  this  independence  the  joint  probability  density  function 
(pdf)  of  x,  y,  u,  v  is  given  by 


p  (x,y,u,v) 

x,y,u,v 


px,u(x-“> 


py,v<>''v) 


(18) 


* 
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where 


P„  (x,u) 


x,u 


XT 


2irabvl-p 


p  (y,v)  = 
y,v 


2ffabVl-p 


(19) 


(20) 


Introduce 


(21) 


Then  (18)  yields 


p  (x,y,u,v) 

Hx,y,u,v 


1 

2  4  2 

(2*)  bV 


— ~-j  (x2  +  y2)  -  ^-(xu  +  yv) 

2b  p  L 


+ 


(22) 


In  order  to  analyze  the  glint  from  the  extended  target  we  require  the 
statistical  properties  of  the  sum  channel  amplitude  and  the  real  and  imaginary 
parts  of  the  complex  monopulse  ratio.  Thus  we  introduce  the  one-to-one 
transformation  of  variables 


5  =  2y-±J~  -  Re  (D/S)  (-“,  ~) 
u  +  v 


(23) 


n 


yu  -  xv 


+  v 


Im(D/S)  (-",  °°) 


(24) 
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C  A  Vu  +  V2  -  |S|  (O.oo) 


9  A  tan-^  ~  ■  phase  of  S  (0,27 r) 


The  pdf  for  the  transformed  variables  may  be  determined  from  the 
transformation  rule 

j-|  j-  g  (Cfn»C»9)  *  px  y  u  v  [*  (C»n » C * 9) » y(5»n* ?»8) ,u(^»n» 5 > 6) » 
v(S,n,?,e>]  J 

where  the  Jacobian  of  the  transformation  Is  given  by  the  determinant 


9x  3x  jbc  _9x 

3£  3n  K  09 


3v  3v  9v  3v 
9n  3£  39 


The  Inverse  transformation  Is  given  by 


=  £(S  cos  9  —  ri  sin  0) 

*  C(S  sin  0  +n  cos  0) 

*  £  cos  0 

*  £  sin  8 


so  that  the  Jacobian  is  explicitly  given  by 

■“  — 
£  cos  9  -X,  sin  0  S  cos  9— r|  sin  0  £(£  sin  0+n  cos  0) 

K  s in  9  £  cos  9  S  sin  0+n  cos  0  £(C  cos  9-n  sin  8) 

0  0  cos  0  -t  sin  0 


cos  0 
sin  0 


-X  sin  0 
£  cos  0 


— ...On..  . . ...v  . 


The  joint  pdf  for  the  transformed  variables  becomes 

,2 


Pr  n  r  e(C*n‘?*6) - 241  e  2b2u2  5  ^  5*1  I  (34) 

(2ir)Vyz  Zb  u  b  * 


9  2  2 

V+rf  -2£a  £+af. 

L  b  b2J 


Since  8  is  uniformly  distributed  in  (0,2ir),  we  have  the  desired  joint 
pdf  for  the  variables  of  interest 


e.n.t 


(S,n,c)  _  £  e 

2irb^y2 


-*[(£)'  -(I)’-] 


(35) 


From  (35)  it  inanediately  follows  that  the  pdf's  for  the  variations  of  £  and  n 
about  their  means  are  identical,  i.e.. 


■  P„(5) 


(36) 


that  each  pdf  is  unlmodal  and  symmetric  about  its  mean,  and  that  these  means 
are 


T  *  P— 
4  b 


n  ■  o 


(37) 

(38) 


Introduce  the  variables 


£  -  Pr 

s  £  -  normalized  Re  (D/S) 

’  y 


(39) 


t  A  2.  normalized  Im  (D/S) 


(40) 


z  A 


yflb 


normalized  S 


(41) 
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which  represent  respectively  the  deviations  of  £  and  n  from  their  means 
normalized  to  y  and  the  sum  channel  amplitude  normalized  to  the  average  sum 
channel  amplitude.  Since  the  Jacobian  of  the  transformation  is  equal  to 

~~  q.  the  joint  pdf  for  the  variables  s,  t,  z  is  given  by 
v^byz 


p.,t,.<8»,:*z) 


3  -z2(s2+t2+l) 
I  z  e 


(42) 


The  variables  s,  t  are  not  independent;  their  joint  pdf  is  given  by 


00 

mJ  T- 


3  e-*V«2«)dz 


T (s2+t2+l) 2 


(43) 


Hence  the  pdf  of  s  is  finally  obtained  as 


Ps(8) 


r  -i— 

J  ir(s2+t2+i)2 


dt 


1  1 


tr 

r-sc 


(s2+l)3/2  11 


2  sec  di  „ 

- 4. 

sec 


2(s2+l)3/Z 


(44) 


“2L 
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The  pdf  of  Re(D/S)  is  given  by  jcf  (39)J  . 

1  1 


2y 

’(5-sf  +1" 

3/2 

v  7  J 

(° s'  ") 


(45) 


and  the  pdf  of  Im(D/S)  is  given  by  £cf  (36)J . 


(46) 
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-3  -3 

Since  the  tells  of  these  pdf's  do  not  decay  faster  than  £  or  n 
second  moments  are  infinite;  as  a  measure  of  spread  about  the  mean 
the  first  absolute  central  moment,  i.e.. 


their 
we  choose 


00 


2 

We  know  that  for  a  Gaussian  pdf  the  variance,  cr_,  is  finite  and  the  first 

u 

absolute  central  moment  is  directly  proportional  to  the  standard  deviation, 
i.e. , 


Thus  the  normalizing  parameter,  y,  defined  by  (21)  corresponds  to  a  standard 
deviation. 

We  next  consider  how  well  thresholding  on  t  and  z  (or  ri  and  £)  reduces 
the  spread  (measured  by  y)  of  the  pdf  of  s  (or  £) . 

The  pdf  of  s  conditioned  on  an  arbitrary  event  A  is 


P 


s 


P(.8ia> _ 

Prob .  (A) 


(49) 


Let  us  now  consider  two  thresholds  T  and  Z  and  a  processor  which  accepts 
measurements  of  Rg(D/S)  only  when  |lm(D/S)|  is  sufficiently  small  and  |s|  is 
sufficiently  large.  In  (49),  the  compound  event  A  is  taken  to  describe  the 
set  of  outcomes 


A  A  {t,z:  | t |  <  T,  z>z) 


(50) 
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t-v 


so  that  (49)  becomes 


N*  (s 


)tj<T,  z>Z) 


✓  1 


w  X 


p(s, t,z)dtdz 


Z  J-  T 


eo  °o  T 

*  f  p ( s, t , z) dtdzds 


/} 


-  H(T,Z) 


(51) 


-T 


Employing  (42)  we  have 
00 

G(s,T,Z)  -  f 

Jz 


2  2 

3  -z  (s  +1) 
z  e 


/"  2  2  , 

2  |  e~z  C  dt 

_  Jo 


dz 


(52) 


In  terms  of  the  error  function  defined  by 


2  f  -t 
erf  x  i  e  dt,  erf(«>) 


(53) 


the  numerator  of  (51)  is  given  by 


G(s 


a 

*T-Z>  ivi 

->7 


2  -z2(s2+l)  ,  N , 

z  e  erf  (Tz)dz 


(54) 


The  denominator,  representing  the  fraction  of  data  which  is  accepted,  is 


H(T, 


00 

Z)  A  j 


G(s,T,Z)ds 


00 

r  -z2 

2ze  erf  (Tz)dz 


■>Z 


(55) 


Integration  by  parts  yields 


H(T,  Z)  -  J"-e~Z  erf  Tzl  +T-2- 

JZ 


-z2(T2+l)  , 
e  dz 


(56) 
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In  terms  of  the  complementary  error  function  defined  by 

at 

2 


erfc  x  A  1-erf  x  ■  — 
■/r 


e-t  dt 


(57) 


The  denominator  of  (51)  becomes 

-Z2  1 

H(T, Z)  -  e  erf(TZ)  + 


erfc 


(Z  >^+1) 


Thus  the  pdf  of  s  conditioned  on  the  event  (t,z: |t|<T,z>z)  is 


(58) 


p.A(s| |t|<T,z>Z)  - 


G(s,T,Z) 

H(T,Z) 


2  f  z2e~z2(s  +1)erf(Tz)dz 

7S  h 


e"Z  erf(TZ)  + 


Vt^+i 


-erfc(Z^) 


(59) 


This  conditional  pdf  has  a  mean  of  zero  and  a  smaller  deviation  about  its 
mean  than  does  the  unconditional  pdf  given  by  (44).  Two  measures  of  this 
deviation  will  be  considered:  the  conditioned  first  absolute  central  moment 
defined  by 


Uc(T,Z)A 


00 

[ 


SIP 


s  A 


(s 


t|<T,z>Z) ds 


AG1(T,Z) 
H(T,  Z) 


(60) 


and  the  conditioned  standard  deviation  defined  by 


<tc(t,z)A 


(61) 
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where 


GX(T,Z) 


g2(t,z) 


|s|g(s,T,Z) ds  - 


s2G(s,T,Z)ds  - 


00 

r 


e-Z  erf(Tz)dz 


00 

r  2 

—  e  z  erf(Tz)dz 
z 


(62) 


(63) 


Recalling  the  relation  (39)  between  s  and  £,  we  note  that  both  the 
conditioned  first  absolute  central  moment  and  the  conditioned  standard  devi¬ 
ation  of  the  pdf  of  C  are  directly  proportional  to  y.  Thus  the  glint  reduction 
in  Re  D/S  accomplished  by  thresholding  will  be  defined  by  the  glint  reduction 
factors 


U  .  H  (T, Z) 

p  (T,Z)  -  G.(T,Z) 

c  J- 


(64) 


.  v. _  4 

^(T.Z, 


1  H  (T,Z) 

2  G2(T,Z) 


h 


(65) 


where  H(T,Z),  G1(T,Z),  G2(T,Z)  are  given  by  (58),  (62),  (63)  respectively. 

We  next  solve  the  following  problems:  choose  the  thresholds  T,Z  so  as 
to  maximize  the  two  glint  reduction  factors  subject  to  a  specified  fraction  of 
data  rejected.  (These  isoperimetric  problems  are  equivalent  to  the  problems: 
choose  the  thresholds  so  as  to  minimize  the  fraction  of  data  rejected  while 
achieving  the  specified  glint  reduction  factors). 

Suppose  we  choose  to  maximize  y"  ^  sub^ect  t0  tbe  constral-nt 


H(T,Z)  -  Hq 


(66) 
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If  we  could  eliminate  T  from  (66)  we  would  have  a  minimization  problem  in  Z 
alone,  l.e.t  we  would  seek  the  Z  which  satisfies 


_d 

dZ 


0 


Unfortunately  an  explicit  expression  for  T(Z)  is  unattainable, 
since  (66)  and  (67)  are  equivalent  respectively  to 


3H  dT  3H  n 
9T  dZ  +  3Z  "  ° 


dT 

we  eliminate  —  and  observe  that  the  extremizlng  thresholds  T,Z 
the  resulting  equation 


3H  ^1  3H  _^fl 
3T  3Z  "  3Z  3T 


Using  (55)  and  (62) 


the  required  partial  derivatives  are  given 


3H  Z  2  -Z2(T2+1) 

9W+1  V? 


3H 

3Z 


-2Ze“Z  erf(TZ) 


-2  -Z 

erf(TZ) 

/w 


(67) 

However , 

(68) 

(69) 

must  satisfy 

(70) 

(71) 

(72) 

(73) 
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3G1  2  1  -Z2(T2+1) 

3T 


When  these  are  employed  in  (70)  we  obtain 


— | erf(TZ)  erfc(Z^+l)-  0 


Each  of  the  two  solutions,  T-°°  and  Z-0,  extremize  the  glint  reduction  factor 
(64)  subject  to  the  constraint  (66).  In  particular  when  we  threshold  only 
on  the  sum  channel  amplitude  (T-00)  we  have  from  (58),  (62) 

-z2 

H(“,Z)  -  e  -  H  (76) 


G^  (“>,  Z)  -  erfc  Z 


Solving  (76)  for  the  threshold,  Z,  in  terms  of  the  fraction  of  data  retained 
we  have 


Z  -  (ln^)1/2 
o 


Thus  one  extremal  of  the  glint  reduction  factor  is  expressed  in  terms  of  the 
fraction  of  data  retained  by 


H  (”,Z)  _ 


U  (»,Z)  G-  (°°,  Z) 

c  l 


[KH 


The  other  extremal  solution  is  obtained  by  thresholding  only  on  the  imaginary 
part  of  D/S,  i.e.,  Z-0.  For  this  situation  (58)  and  (62)  yield  respectively 


H(T,0) 


/t^+i 


r 


Gi  (T-°)  -  x 


e~Z  erf(Tz)dz 


J0 


(81) 


The  Integral  may  be  evaluated  by  noting  that 


dG 

ar<T-0) 


? 


2^ 

¥ 


2.e-l2<l2+l>d2  -  I  1 


-'O 


l^+l 


(82) 


Since 


G1(0,0)  -  0 


(83) 


we  have 


G1(T,0) 


(84) 


From  (80)  we  find  the  threshold  in  terms  of  the  fraction  of  data  retained 

T  -  -j£==  (85) 

/l-Ho2 

hence  the  second  extremal  of  the  glint  reduction  factor  expressed  in  terms  of 
the  fraction  of  data  retained  is  given  by 


VT fi) 


H  .  H  (T»0)  _ 
G- (T,0) 


Ho 

(86) 

to) 
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For  any  other  combinations  of  threaholds  which  satisfy  the  constraint 


-Z2 


erf(TZ)  + 


erfc 


■  Ho 


(87) 


we  obtain  the  glint  reduction  factor 


V  m  H  (T,  Z) 

U  (T,Z)  G.(T,Z) 

c  l 


Ho 


_  2 

e-Z  erf(Tz)dz 


(88) 


which  is  Intermediate  between  the  extremal  values  given  by  (79)  and  (86). 

The  second  measure  of  glint  reduction  is  given  by  (65).  Using  (63) 
we  find  the  partial  derivatives. 


8G2  -1  -Z2 

-gf  -  Z  erf(TZ)  (89) 


3G2  -  — =erfc(Z*45+l)  (90) 

9T  /r2+l 


The  equation  for  the  extremizing  thresholds, 


9H  ^2  9H  ^2 
9T  9Z  ~  9Z  9T 


(91) 


thus  yields 


.Z2 

e  erf(TZ) 
ZCT2-*-!)^2 


|  [l-2Z2(T2+l)] 


erfc(Z«42+l)  +  Zv'P+l  e_Z  (T  +1) 


0 


Clearly  one  solution  Is  given  by  T-”,  i.e.,  threshold  on  only  the  sum  channel 
amplitude.  For  this  situation  (58)  becomes 


H(»,Z) 


(93) 


and  (63)  becomes 


(94) 


where  E^(.)  is  the  exponential  Integral.  Note  that  this  becomes  unbounded  as 
Z-*C.  This  is  to  be  expected  since  the  situation  T"00,  Z«0  corresponds  to  Ho*l, 
i.e.,  to  a  processor  which  accepts  all  measurements  of  Re  D/S.  The  pdf  of  its 
output  is  given  by  (45)  which  as  we  have  already  noted  has  infinite  variance. 
From  (93) 

Z  -  (ln^)1/2  (95) 

so  that  an  extremal  of  the  second  glint  reduction  factor  is 


(96) 


We  now  prove  that  this  is  the  only  extremal.  Note  that  any  other  combination 
of  thresholds  which  provide  an  extremal  must  satisfy  [cf.(92)]. 


(l-2x2)  erfc(x)  +  —  xe  x  -  0 

/F 


(97) 


where 


x 


>  0 


(98) 
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Clearly  x  ■  is  not  a  solution.  Hence  we  seek  values  of  x  for  which 

2 


J  2  -x 

e  dt  +  -SI—-  -  0 
1-2  x 


(99) 


Since 


f  (o)  -  —  >  0  (100) 

and 

,  2  -x2 

f'(x)  -  ■-■*-= . -  >  0  (101) 

(l-2xV 

we  conclude  that  f(x)  never  vanishes  for  positive  x  and  hence  that  (97)  has 
no  positive  roots. 

Thus  the  only  extremal  of  —5—^ -  subject  to  the  constraint  (66)  is 

/k(T-2> 

given  by  (96).  When  we  threshold  on  T  alone  we  do  not  obtain  an  extremal, 
however  the  glint  reduction  factor  for  this  case  is  obtained  by  noting  that 
(58)  and  (63)  yield  respectively 


H(T,0) 


(102) 


g2(t,o) 


J  ^e  z  erf(Tz)dz 
0 


To  evaluate  the  integral  we  note  that 


(103) 


G2(0,0)  -  0 


(104) 
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00 


(105) 


hence 


g2(t,o> 


tan-1T 


0 


(106) 


Note  that  G2(T,0)  and  hence  <Jc(T,  0)  has  a  logarithmic  singularity  at  T«°°. 
Using  (102)  this  is  given  in  terms  of  the  fraction  of  data  retained  as 


g2(t,o) 


(107) 


so  that  finally  the  (non-extremal)  glint  reduction  factor  obtained  by 
thresholding  on  T  alone  is  given  by 


V 

t r  H  (T, 0) 

1/2 

TTHO 

/|a  (T,0) 

TT  C 

2  G2(T,0) 

1  (  l4fl0\ 
L1"  (,l-Ho)J 

(108) 


III.  DISCUSSION  OF  RESULTS 

The  results  of  this  analysis  are  shown  numerically  in  the  accompanying 
graphs.  Figures  3  and  4  show  respectively  graphs  of  the  first  glint  reduction 
factor  —  ^  and  the  fraction  of  the  data  rejected  1-H(T,Z)  vs  sum  channel 

amplitude  threshold  Z  for  values  of  fixed  thresholds  on  |lm^[  equal  to  .05, 

.1,  .3,  .5,  1,  1.5,  2,  1000.  Figure  5  shows  the  cross  plots  of  - — .  and 

y^\ l) 

includes  the  two  extremals.  It  is  seen  that  for  a  given  fraction  of  rejected 
data,  maximum  glint  reduction  is  achieved  by  thresholding  on  Z  alone.  Along 
each  of  the  intermediate  curves  T  is  held  constant  and  Z(T)  is  chosen  to  yield 
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the  specified  value  of  1-Ho.  Figure  6  shows  two  examples  of  the  second  glint 
reduction  factor  vs  1-Ho.  The  upper  curve  is  the  sole  extremal;  the  lower 
curve  is  a  special  threshold  selection  (z*0)  which  corresponds  to  thresholding 
only  on  |lm(D/S)|.  As  for  the  first  glint  reduction  factor,  thresholding  on 
Z  alone  provides  better  performance. 

-1/2 

The  dashed  curve  in  Figure  6  is  a  plot  of  the  function  H  vs  1-H. 

It  pertains  to  independent  samples  from  any  stationary  random  process  with 
finite  variance.  It  shows  the  reduction  in  the  standard  deviation  of  the 
sample  mean  if  all  samples  are  retained.  This  potential  reduction  is  never 
realized  when  the  process  is  Re  D/S  since  the  unthresholded  samples  do  not 
have  finite  variance.  Nevertheless  we  conclude  that  for  almost  the  entire 
domain  of  the  graph,  we  achieve  a  greater  measure  of  glint  reduction  by 
retaining  only  the  fraction  H  of  the  measurements  whose  sum  channel  amplitude 
exceeds  the  threshold  Z  than  by  forming  a  sample  mean  from  all  the  Independent 
measurements. 
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1— H(T,Z) 


Fi*ir*  4.  Fraction  of  data  rejected  vt.  eum  channel  th refold. 
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Figure  6.  Glint  reduction  factor  vi.  fraction  of  data  refactad. 
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Figure  6.  Glint  reduction 
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The  preceding  analysis  gives  the  glint  reduction  when  thresholds  are 
applied.  In  order  to  find  the  actual  glint  before  thresholding,  the  elements 

of  the  covariance  matrix  must  be  determined.  These  are  given  by  (15),  (16), 

—  IT 

and  (17).  Assume  that  |A(iJ>)  \  is  independent  of  0,  i.e.,  assume  a  uniform 
distribution  of  scatterers  across  the  target.  Then 


(A-2) 


2pab 


( 


v! 


■  5 
K-J 


(A-3) 


We  make  the  observation  that  the  integral  representation  of  the  Bessel 
function  of  the  first  kind  of  integer  order 

J„(.>  *  f  e‘<Z  *♦ 

J—TT 


is  the  n^  Fourier  coefficient  in  the  expansion  of  eiz  sln^ 
the  Anjer-Weber  representation 


.2rrd 
e1  X 


sinijj 


Thus  we  obtain 


(A-5) 
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Integrating  over  the  target  extent  we  have 


rVf 


C  o  2  t  ■  i  2 

J  { .*¥  J»(2f )  .‘“♦o  j  4  m**i 


(A-6) 


T^»e  elements  of  the  covariance  matrix  become 


*4*  K(¥) +  2  J  j2„  (¥)“«"  ♦.  ¥¥] +  |ns|2  <A'7) 


■  K(¥)  - 2  t  J2»  (¥)“'2"*  ¥¥] + 1": 


(A-8) 


2pab 


Bin(2n+l)4»o 


sln(2n4l)6/2  1 

(2n+l)6/2  J 


(A-9) 


where  we  have  employed 

J_n^z)  ”  (-I)”  Jn(z)  (A-10) 

Thus  for  specified  antenna  dimension,  RF  wavelength,  target  extent  and  average 
noise  powers,  the  spread  of  the  pdf  of  Re  D/S  before  thresholding  may  be 
determined. 


DISTRIBUTION 


No.  of 
Copies 


IIT  Research  Institute  1 

ATTN:  GACIAC 
10  Uest  35th  Street 
Chicago,  IL  60616 

DRSMI-R,  Dr.  McCorkle  1 

-RE  1 

-RE,  Mr.  Pittman  1 

-REL  2 

-REG  1 

-REM  1 

-RER  20 

-RPT  (Record  Copy)  1 

-RPR  15 

DRSMI-LP,  Mr.  Voigt  1 


31 


